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We use scanning tunneling microscopy to visualize and thermal desorption spectroscopy to quan¬ 
titatively measure that the binding of naphthalene molecules to graphene (Gr), a case of pure van 
der Waals (vdW) interaction, strengthens with n- and weakens with p-doping of Gr. Density func¬ 
tional theory calculations that include the vdW interaction in a seamless, ab initio way accurately 
reproduce the observed trend in binding energies. Based on a model calculation, we propose that 
the vdW interaction is modified by changing the spatial extent of Gr’s tt orbitals via doping. 


One of the key properties of graphene (Gr) is the 
wide-range tnnability of its Fermi level and correspond¬ 
ing charge carrier concentration, either by a gate elec¬ 
trode [1], substitutional doping [2], adsorption [3, 4], or 
charge transfer from a supporting material or intercala¬ 
tion layer [5-8]. The tnnability of the Fermi level through 
the otherwise rigid band structure results from the mate¬ 
rial being atomically thin and having a negligible density 
of states near the Dirac point. 

In recent years, interest has arisen in using this tunabil- 
ity to control adsorption: For the case of ionic adsorbates, 
Brar et al. [9] demonstrated a dependence of the ioniza¬ 
tion state of a Go adatom on the Gr Fermi level position, 
and Schumacher et al [10] found a doping-dependent 
binding energy Ejj of ionic adsorbates to Gr, with a shift 
in El) on the order of the shift in the Fermi level in¬ 
duced by doping. For the case of radicals and based on 
ab initio calculations, Wehling et al. [11] predict doping- 
dependent adsorbate phase transitions for hydrogenated 
as well as fluorinated Gr, while Huang et al. [12] find a 
stronger binding of isolated H radicals for larger magni¬ 
tudes of doping. 

For the case of van der Waals (vdW) interaction, the 
effect of the Gr doping level on the binding energy of 
adsorbates has not yet been explored. This is surpris¬ 
ing, given that the adsorption of simple hydrocarbons to 
graphite or Gr has been used as a model system to study 
vdW interactions [13-16]. Here, we investigate this case 
with the help of epitaxial Gr on Ir(lll), which can be 
doped from the backside by intercalation of highly elec¬ 
tropositive (e. g. Gs, Eu) or -negative (e.g. O) ele¬ 
ments into its interface with the substrate, while Gr’s 
other side remains available for the adsorption experi¬ 
ment itself. This strategy not only enables us to achieve 
large Fermi level shifts on the order of ±leV, but also 
to visualize doping-induced binding energy differences by 


making using of intercalation patterns [10]. Naphthalene 
is chosen as a test molecule, since its binding to Gr is a 
pure vdW case studied previously, both experimentally 
[14] and theoretically [13]. Our experiments are comple¬ 
mented by density functional theory (DFT) calculations 
that include the vdW interaction in a seamless, ab initio 
way (for a recent review, see Ref. [17]). 

For this paradigmatic case we find in excellent agree¬ 
ment of experiment and theory an increase of the vdW 
binding energy when changing from p- to n-doping. We 
identify the mechanism of this vdW binding energy ma¬ 
nipulation, which is likely to hold for a broad class of 
molecular adsorbates and a variety of 2D materials. 

The experiments were performed in two variable- 
temperature scanning tunneling microscopy (STM) sys¬ 
tems with base pressures below 1 x 10“^^mbar. Ir(lll) 
was prepared by cycles of sputter-anneal, yielding clean 
terraces with sizes on the order of 100 nm. Well ori¬ 
ented Gr was prepared by room temperature ethylene ad¬ 
sorption till saturation, thermal decomposition at 1500 K 
and subsequent exposure to 5 x 10“^ mbar of ethylene at 
1170K [18]. Using 600s ethylene exposure time, a per¬ 
fectly closed, 1 monolayer (ML) Gr layer was achieved, 
while with 180 s a partial ^ 0.9 ML Gr layer resulted. 
For n-doping, Eu was intercalated under 1 ML Gr on 
Ir(lll) at 550K by exposure to high purity Eu [19] va¬ 
por. Eor p-doping, oxygen was intercalated till satura¬ 
tion under 0.9 ML Gr on Ir(lll) at 530K using a lo¬ 
cal oxygen pressure of about 2.5 x 10“^ mbar for 600 s. 
Naphthalene exposure [20] was controlled by a variable 
leak valve. STM images were taken at 35 K and digi¬ 
tally post-processed with the WSxM software [21]. Ther¬ 
mal desorption spectroscopy (TDS) was performed with 
a sample heating rate of 1 Ks“^ by monitoring desorb¬ 
ing naphthalene (molecular mass m = 128 u). The sam¬ 
ple temperature measurement during TDS with a K-type 
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thermocouple has an absolute error below 3 K with a re¬ 
producibility better than 1 K. 

The DFT calculations were carried out in the gener¬ 
alized gradient approximation (GGA) [22] with a kinetic 
energy cut-off of 500 eV using the projector augmented 
wave method (PAW) [23] as implemented in the VASP 
code [24, 25]. The ground-state geometry and the cor¬ 
responding electronic structure were investigated using 
the non-local correlation energy functional vdW-DF2 [26] 
with a revised [27] Becke (B86b) exchange energy func¬ 
tional [28]. The unit cell comprised (10x10) graphene 
unit cells on 4 layers of (9x9) Ir(lll), either pristine or 
(Eu/0)-intercalated. 

We start with a discussion of the STM preparations 
and measurements. Upon intercalation of sub mo nolayer 
amounts of Eu, a complex pattern of stripes and islands 
of Eu-intercalated Gr (high/bright) in coexistence with 
non-intercalated Gr (low/dark) emerges as shown in the 
STM topograph in Eig. 1(a) and discussed in Ref. [29]. 
Within the intercalated patches, Eu forms a p(2 x 2) 
superstructure and leads to local doping of graphene by 
— 1.36 eV, while the surrounding pristine Gr areas remain 
hardly doped at +0.1 eV [10]. 

After deposition of submonolayer amounts of naph¬ 
thalene at T < 50K and subsequent brief annealing to 
T = 150 K, the sample is cooled down again for STM 



FIG. 1. (color online) (a-c): STM topographs imaged at 
35K. (a): Gr/Ir(lll) after intercalation of about half 
a monolayer coverage of Eu ((39nm)^). (b): Like (a), 

but with additional submonolayer coverage of naphthalene 
((90 nm)^). (c): Molecular resolution of naphthalene island 
((18 X 11) nm^). (d): Schematic cross section corresponding 
to profile in (e). (e): Height profile along the white line indi¬ 
cated in (b). 



temperature T [K] 


FIG. 2. (color online) Mass spectrometer signal of desorbing 
naphthalene (m = 128 u) while heating the sample with a 
ramp rate of IK/s in dependence of intercalant for the cov¬ 
erages of 1 ML and ~ (2.3 ± 0.2) ML each. 


imaging. The corresponding Fig. 1 (b) displays addi¬ 
tional naphthalene islands (cyan encircled) only on top 
of Eu-intercalated patches. The profile in Eig. 1 (e) taken 
along the white line in Eig. 1 (b) shows that the appar¬ 
ent height of the Eu intercalation island itself stays un¬ 
changed at 2 A, while the naphthalene islands have an 
apparent height of about 1.3... 1.5 A on top of the Eu- 
intercalated Gr. The schematic cross section of Eig. 1 (d) 
visualizes the morphology underlying the height pro¬ 
file. Higher-resolution imaging displayed in Eig. 1 (c) 
reveals molecular resolution of the naphthalene adsor¬ 
bate layer. Our interpretation is as follows: During an¬ 
nealing, the molecules are highly mobile, diffuse over the 
surface while sensing local adsorption energy differences. 
As the sample is cooled down again, the arrangement is 
frozen and the preferential coverage of intercalated areas 
by molecules then directly reflects the preferential bind¬ 
ing to these areas. 

To strengthen our case for the generality of this result, 
we repeated the experiment exchanging naphthalene with 
benzene, and found qualitatively the same results (com¬ 
pare [30]). 

In order to quantify the difference in E^, observed in 
STM, we investigated thermal desorption of naphthalene 
from homogeneous Gr layers, z.e., which are either en¬ 
tirely pristine or entirely Eu p(2 x 2)-intercalated. 

Eigure 2 shows the desorption traces for naphtha¬ 
lene coverages of 1 ML and (2.3 ± 0.2) ML. The low- 
temperature peak corresponds to desorption from the 
naphthalene multilayer and is zero-order as seen in the 
exponential leading edge followed by a sharp drop in in¬ 
tensity, while the high temperature peak is almost sym¬ 
metric, indicative of first-order desorption. Eor pristine 
graphene, the maximum of the monolayer peak is at 
Tmax = 236 K. These findings agree in the limits of 
error with those of Zacharias et al.[14] for naphthalene 
desorption from graphite, a hint that the influence of the 
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Ir substrate on the binding of naphthalene to Gr is neg¬ 
ligible. 

Upon intercalation of Eu (2x2) (associated with a 
— 1.36eV Gr doping level), the monolayer peak is shifted 
up to 244 K, evidence for increased while the po¬ 
sition of the multilayer peak is the same as seen in its 
intercept with the desorption rate at the maximum of 
the monolayer peak (horizontal dashed lines) at 196 K 
(left vertical dashed line). As the multilayer peak is not 
expected to shift due to the diminishing effect of the sub¬ 
strate with increasing molecular film thickness, its pres¬ 
ence provides an integrated temperature calibration. 

As shown in Fig. 2, by TDS we also investigated O- 
intercalated Gr, for which the monolayer Tj^ax is shifted 
down to 233 K. This case would have been difficult to 
analyze by STM, as the intercalation is not regularly 
patterned as for the case of Eu. TDS of naphthalene 
adsorbed to two other Gr intercalation systems [Cs, Eu- 
(\/3 X \/3)] are analyzed in the Supplemental Material 
[30]. Altogether, we find a strictly monotonic increase of 
\Ei)\ with the Fermi level moving up in the bandstructure 
from p- to n-doping. 

To extract E^, from the desorption temperatures, we 
employed the Redhead method using the previously de¬ 
termined value of the pre-exponential factor z/ = 5 x 
1016±2 hz [ 14 ]. The left part of Table I lists the exper¬ 
imentally determined E^^"^ alongside the doping levels 
E^"^, which are given in terms of the position of the Dirac 
cone relative to the Fermi level as determined in angle- 
resolved photoemission (ARPES) measurements. The 
experimental error in is about 100 meV, and domi¬ 
nated by the error of the pre-exponential factor. How¬ 
ever, as the pre-exponential factor acts equally in the 
calculation of all the errors of the differences in E^, 
for the differently doped Gr are dominated by the re¬ 
producibility of the temperature measurement, and are 
therefore much smaller, only about 4meV. 

In order to better understand the origin of the ob¬ 
served change in Ef) upon intercalation, we have con¬ 
ducted DFT calculations of the naphtalene molecule on 
non-intercalated as well as O- and Eu-intercalated Gr. 
Figure 3 (a) and (b) show the geometry of the calcula- 


experiment theory 



^exp rji 

D -^max 

[eV] [K] 

^exp 

[eV] 

TTitot 

ipDFT ipNL 

[eV] 

TTitheor 

Gr/O/Ir 

+ 0 . 68 ® 233 

- 0.808 

- 0.806 

0.366 - 1.172 

+ 0.88 

Gr/Ir 

+ 0 . 10 '^ 236 

- 0.819 

- 0.852 

0.556 - 1.408 

+ 0.15 

Gr/Eu/Ir 

- 1 . 36 '' 244 

- 0.848 

- 0.865 

0.736 - 1.601 

- 1.21 


^Ref. [31] (also Refs. [7] and [32]) ^ Ref. [33] Ref. [8] 


TABLE I. Doping levels Ed, desorption peak maximum Tmax 
and binding energies Eb from experiment and theory. 
is given assuming jy = 5 x 10^® Hz (without error as it acts 
equally on all values). 



FIG. 3. (color online) Top view of relaxed adsorption 
geometries of naphthalene on (a) Gr/0/Ir(lll) and (b) 
Gr/Eu/Ir(lll). In (c) and (d), the corresponding plots of the 
charge density difference for each system are also presented. 
Color scale range is ±5 x 10“^ e/ao with ao the Bohr radius. 

tions for the intercalated systems, while the right part 
of Table I provides the calculated values for doping level 
and binding energy 

The calculated doping level has been determined 

from the shift of Gr’s Dirac point against the Fermi level 
and should thus be directly comparable to the value ex¬ 
tracted from ARPES. Indeed, we find that the calculated 
doping levels are in reasonable agreement with the 

experimentally determined in particular confirming 

n-doping of Gr by Eu and p-doping of Gr by O intercala¬ 
tion. Qualitatively, the doping is also seen in the charge 
density difference plots in Fig. 3 (c) and (d), which give 
a direct image of charge accumulation (for Eu) and de¬ 
pletion (for O) on the Gr site. These plots also visualize 
that charge transfer between molecule and Gr is absent 
(a Bader analysis [34] yields a charge transfer of less than 
0.03 electrons), which is a precondition for naphthalene 
adsorbed to Gr being a model for vdW interaction even 
when Gr is supported on a substrate and doped by inter¬ 
calation. 

Next, we look at the theoretical binding energy, de- 
fined as El°^ = El°l - {El°l{ + where El°l is 

the total energy of the molecule-surface system, E^^^^ de¬ 
notes the total energy of the surface (i.e., the Gr/Ir(lll) 
and Gr/Eu/Ir(lll) systems) and E^^oiec corresponds to 
the total energy of the naphthalene molecule in the gas 
phase. As obvious from Table I, the calculations repro¬ 
duce the experimentally observed increase in \Eb\ rather 
well when going from Gr/0/Ir(lll) via Gr/Ir(lll) to 
naphthalene/Gr/Eu/Ir(lll). Yet, considering absolute 
numbers, two points have to be noted: Firstly, while Gr 
is essentially freestanding when intercalated with O or 
Eu [7, 8], Gr on bare Ir(lll) is slightly hybridized with 
its substrate in specific areas of its moire [35]. This in¬ 
duces a certain modulation of E^ and thus renders the 
non-intercalated case less reliable. Secondly, the excel¬ 
lent match of the absolute values (within < 25 meV) of 
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FIG. 4. (color online) (a) Schematic band structure near 
the Dirac point and representation of the tt orbital exten¬ 
sion for p-, undoped and n-doped graphene, (b) In-plane 
a:p-integrated 7r-state charge density for (10x10) unit cells 
of freestanding Gr as a function of ^-direction, (c) Zoom into 
the low-density tail of (b). 

the binding energies must be considered fortuitous, as Ei^ 
extracted from TDS is expected to carry a larger system¬ 
atic error. Unaffected by these two caveats, however, is 
the difference in Eij between the O- and Eu-intercalated 
cases, which matches reasonably well being ^ 60 meV in 
the calculation and ~ 40meV in the experiment. 

To gain an insight into the role of the vdW interactions 
in these systems, the total binding energy E\f^ was de¬ 
composed into a DFT contribution E^^^ and a non-local 
(NL) one Ef^ (for details, see Supplemental Material 
[30]) as seen in Table 1. The data show that for all systems 
considered in our study, the DFT contribution E^^^ to 
is positive while the NL one E^^ is negative indi¬ 
cating that the vdW interactions are indeed the driving 
force leading to a stable molecular adsorption. Further¬ 
more, the magnitude of the attractive NL contribution 
Ef'^ to the naphthalene binding significantly increases 
in the series from the Gr/0/Ir(lll) (1.172 eV) via the 
Gr/Ir(lll) (1.408eV) to the Gr/Eu/Ir(lll) (1.601 eV) 
system. 

To find out the reason for the modification of the vdW 
binding upon intercalation, we have conducted calcula¬ 
tions on a freestanding Gr sheet of (10 x 10) unit cells. 
It was doped by directly adding or removing 0.05 elec¬ 
trons per C atom. To keep the calculated unit cell on the 
whole charge-neutral, a balancing background charge had 
to be introduced. Eigure 4 (b) shows the resulting xy- 
integrated charge density for all tt orbitals corresponding 
to the differently doped cases. As qualitatively illustrated 
in Eig. 4 (c), the p- or n-doping of Gr results in a tt charge 
density distribution that is less (p-doped Gr) or more (n- 
doped Gr) spatially extended as compared to that of the 


undoped Gr. A less (more) spatially extended tt charge 
density corresponds to less (more) polarizable C atoms in 
the Gr layer and therefore will govern a weaker (stronger) 
vdW binding of the tt adsorbates. This picture is derived 
from the behavior observed for the vdW bonded noble- 
gas dimers from He to Kr [36, 37]. It is backed up quan¬ 
titatively also by a recent analysis of Borland et al. [38], 
who found for benzene on Gr, that most of the non-local 
energy in the vdW-DE-type functional arises from the 
extended regions of low charge density. 

The increase in the non-local contribution with elec¬ 
tron doping in Table I is mostly counteracted by the si¬ 
multaneous increase in the repulsive DET contribution, 
such that the total binding energy Ef^^ of the molecule to 
the surface is modified only slightly. Generally in a pure 
vdW system, the attractive vdW forces have to be bal¬ 
anced by a counteracting Pauli repulsion force in any sta¬ 
ble adsorption geometry. If one approximates the vdW 
contribution by the non-local correlation, one might ten¬ 
tatively attribute the increase in the repulsive DET term 
to an additional Pauli repulsion as Gr’s charge density 
becomes more extended due to n-doping. 

In view of the generality of the mechanism proposed 
above, we speculate that also substitutional doping or 
tuning of the Gr Eermi level through a gate electrode 
will affect the vdW contribution of adsorbate binding in 
a similar way. 

In conclusion, in excellent agreement of experiment 
and theory we find an increase in the vdW binding en¬ 
ergy of naphthalene to Gr by about 5%, when chang¬ 
ing from p- to n-doping in the experimentally accessible 
range. Based on our DET calculations, the Gr tt orbitals 
become spatially more extended with n-doping, making 
the space of low electron density between Gr and the 
molecule more polarizable, which in turn gives rise to a 
stronger vdW binding. The effect is robust with respect 
to changing either the intercalant or the molecule and 
must therefore be assumed to provide a broadly applica¬ 
ble strategy to manipulate van der Waals binding in tt-tt 
systems. 
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This Supplemental Material presents thermal desorp¬ 
tion spectroscopy of naphthalene on additional interca¬ 
lation systems, scanning tunneling microscopy (STM) of 
benzene on intercalation-patterned graphene (Gr), and 
a description of the binding energy evaluation from our 
density functional theory (DFT) calculations. 


FULL SET OF DESORPTION TRACES 
FOR MORE INTERCALATION STRUCTURES 


In the paper, we presented desorption of naphthalene 
from pristine Gr/Ir(lll) and intercalated either with 
oxygen to saturation or Eu in a p(2 x 2) superstructure, 
and only two selected coverages of naphtalene for reasons 
of clarity. 

However, we also investigated two other Gr intercala¬ 
tion structures: (1) Gs adsorbed in a (a/3x v^)R30° with 
respect to Ir, prepared by room temperature deposition 
of Gs from commercial dispensers [1] till saturation, and 
(2) Eu adsorbed in a (\/3 x \/3)R30° with respect to Gr, 
prepared by deposition of Eu at 720 K until saturation. 

The full set of all our measured desorption traces for 
the in total five different intercalation systems and many 
different coverages each are shown in Eig. SI. The cover¬ 
ages were determined by integrating the area under the 
background-substracted spectra. 

As all intercalation structures are associated with a 
different doping level (see SI for references), we obtain 
one point for each in our plot of desorption temperature 
vs. doping level in Eig. S2. While our expectation of 
a monotonous increase of naphthalene’s binding energy 
with increasing electron transfer into graphene is con¬ 
firmed, the effect seems nonlinear and a simple functional 
dependence cannot be inferred from our present data. 



FIG. SI. Mass spectrometer signal of desorbing naphthalene 
{m — 128 u) while heating the sample with a ramp rate of 
1 K/s for different intercalants with naphthalene coverages in 
the range up to a few monolayer (ML). 
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Ed 


Tmax 

^exp 


[eV] 


[K] 

[eV] 

Graphite 

0 


235 

[2] 0.816 

Gr/Ir 

+0.10 

[3] 

236 

0.819 

Gr/Eu p(2 X 

2)/Ir -1.36 

[4] 

244 

0.848 

Gr/Eu {\/3 X 

V^)R307Ir -1.43 

[4] 

246 

0.848 

Gr/Cs {Vs X 

V3)R307lr -1.13 

[5] 

239 

0.848 

Gr/O/Ir 

+0.68 

[6], [7, 8] 

233 

0.808 


TABLE SI. Experimental doping levels Ed, desorption peak 
maxima and resulting binding energies Eh for desorption of 
naphthalene from differently doped Gr. is given assum¬ 

ing z/ = 5 x 10^® Hz [2]. 



Doping level Eq [eV] 

EIG. S2. Desorption temperature of naphthalene vs. 
graphene doping. Each data point corresponds to a certain 
intercalation structure. 


BENZENE ADSORPTION ON 
INTERCALATION-PATTERNED GR 


In the paper, we described that on intercalation- 
patterned Gr, naphthalene islands are preferentially lo¬ 
cated on the En-intercalated areas. We have repeated 
the experiment with benzene as a test of generality, and 
found the same result, as shown in Fig. S3 (a)-(d). 

However, in the STM topograph in Fig. S3 (a), another 
peculiarity is visible: It stands out that benzene islands 
fill an intercalation island always either completely or 
not at all. Furthermore, larger intercalation islands are 
predominantly covered with benzene islands (the three 
largest islands in the image are all covered), while smaller 
intercalation islands are less often and thin stripes never 
covered. Our interpretation is that this is a confined 
Ostwald ripening process: As larger islands have smaller 
edge-to-area ratios, they are energetically more favor¬ 
able. During the annealing, thermal energy is sufficient 
to make the molecules mobile also over non-intercalated 
ares. Then, larger benzene islands grow at the expense 
of smaller ones until the area of preferential binding on 
the intercalation island is completely filled. 



0 30 60 90 120 



FIG. S3, (a) and (d): STM topographs imaged at 35K. (a): 

Partially Eu-intercalated Gr/Ir(lll) with submonolayer cov¬ 
erage of benzene ((120 nm)^). (b): Schematic cross section 
corresponding to the line profile in (c). (c): Height profile 
along the white line indicated in (a), (d): Molecular resolu¬ 
tion of a benzene island ((20 x 12) nm^). 


EVALUATION OF THE BINDING ENERGY 
FROM THE DFT CALCULATIONS 

In DFT [9], the total energy of the physical system 
under consideration is a functional of the electron charge 
density n(r) and can be expressed in the Kohn-Sham 
formalism [10] as 

E[nY'^^ = Ts[n] + F’extM + Eh[ti] + Ex[n] Ec[n]. (1) 

where Ts[n] describes the kinetic energy, F^ext[^] is the 
potential energy due to the external potential (z.e., gener- 
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ated by nuclei, etc.), Enin] corresponds to the potential 
energy due to the repulsive Coulomb electron-electron 
interaction and Ex[n] and Ec[n] represent the exchange 
and correlation energy, respectively. In particular, in the 
commonly used PBE exchange-correlation functional [11] 
the correlation part is given by a sum of a LDA and a 
semi-local term 

EM-E];^^[n]+Ef[nl ( 2 ) 

while in the vdW-DF formalism [12] the semi-local com¬ 
ponent E^^[n] of the correlation energy is replaced by a 
non-local term £’^^[n] 

EM=E^^^[n] + Ef^[n]. (3) 

In consequence, the total energy of the system can be 
cast into the following expression 

E [n] = E[nf^'^ + Ef^[n], (4) 

where 

E[nf^^ = T,[n] + E,^t[n] + Enin] + E,[n] + E^^^[n]. 

( 5 ) 

Finally, we employed Eq. 5 together with the definition of 
the binding energy of the molecule-surface system 

El°^ = ^molec+surf “ (^molec + ^suJf) tO decompose the 
El°^ into a DFT and a NL part 

Etot ^ ^DFT ^ (g) 
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